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ABSTRACT 

In this work, review of various well-known 
algorithms for image denoising is carried out & their 
performances with their methodologies are 
comparatively assessed. A new algorithm based on 
the orthonormal wavelet transform (OWT) is 
developed. In this work images corrupted by AWGN 
are denoised. Simulation results shows that proposed 
method using Orthonormal wavelets for different 
values of noise Standard Deviation a in dB 

outperforms other available methods. Also Coiflet 
Wavelet performs better than Symlet, Haar & 
Daubechies wavelets. The proposed Orthonormal 
wavelet transform (OWT) method has minimum 
Mean Square & highest PSNR with Coif let wavelets. 
Simulation results shows that denoised image is 98.29 
% similar for 5 dB noise standard deviation and 
84.42% similar for 30 dB noise standard deviation. 
The proposed method has 1.35 dB & 4% 

improvement for a =10 dB, 2.08 dB & 7% 

improvement for a =20 dB & 2.26 dB & 9% 

improvement for a =30 dB as compared to denoising 
with two thresholds for edge detection [1], 

Keyword: A WGN, Image Denoising, Noise, 

Thresholding, DWT, OWT, PSNR, MSE, SSIM. 

I. INTRODUCTION 

Image denoising is often used in the field of 

photography or publishing where an image was 
somehow degraded but needed to be improved before 
it can be printed. For this type of application, we need 
to know something about the degradation process in 
order to develop a model for it. When we have a 
model for the degradation process, the inverse process 
can be applied to the image to restore it back to the 
original form. Denoising is an important part of many 
image processing algorithms [1]. It is also one of the 


main challenges in image processing. Noise reduction 
improves perception of images and usually results in 
better performance for different image analysis and 
computer vision algorithms. In our model we suppose 
that an image is distorted by additive white (i.e. 
uncorrelated) noise. One-level discrete wavelet 
transform (DWT) might be treated as a separation of a 
signal into higher and lower frequency bands with 
further down-sampling components [1], It is possible 
to apply such discrete one-level wavelet transform to 
lower frequency component recursively to perform 
two-level transform. If applied K times the frequency 
separation procedure results in K-level DWT [1], 

Wavelet transform describes both frequency (like 
Fourier transform) and time characteristics of signal. 
It allows finding edges directly in the wavelet 
transform domain. Two-dimensional discrete wavelet 
transform is implemented by ID transform applied to 
rows and columns of 2D-signal (in our consideration, 
a grayscale image) separately and subsequently. As a 
result of one-level 2D-DWT we get the wavelet 
spectrum with four squared discrete sets called sub 
bands: one with lower frequency components and 
three sub bands with higher frequency components 
associated with horizontal, vertical and diagonals 
edge directions. 


Figure 1: Original image Discrete Wavelet Transform 
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II. Discrete Wavelet Transform (DWT) 

Wavelets are mathematical functions that analyze data 
according to scale or Resolution. They aid in studying 
a signal in different windows or at different 
resolutions. For instance, if the signal is viewed in a 
large window, gross features can be noticed, but if 
viewed in a small window, only small features can be 
noticed. Wavelets provide some advantages over 
Fourier transforms. For example, they do a good job 
in approximating signals with sharp spikes or signals 
having discontinuities. Wavelets can also model 
speech, music, video and non-stationary stochastic 
signals. Wavelets can be used in applications such as 
image compression, turbulence, human vision, radar, 
earthquake prediction, etc. The term “wavelets” is 
used to refer to a set of Ortho-normal basis functions 
generated by dilation and translation of scaling 
function cp and a mother wavelet \|/. The finite scale 
multi resolution representation of a discrete function 
can be called as a discrete wavelet transform. DWT is 
a fast linear operation on a data vector, whose length 
is an integer power of 2. This transform is invertible 
and orthogonal, where the inverse transform 
expressed as a matrix is the transpose of the transform 
matrix. The wavelet basis or function, unlike sine and 
cosines as in Fourier transform, is quite localized in 
space. But similar to sine and cosines, individual 
wavelet functions are localized in frequency. The 
Ortho-normal basis or wavelet basis is defined as: 

<PiJ4d M H 2i<p(¥x - fc) 

The scaling function is given as: 

a (a) = 2?<p(2ix — &} 

Where ip is called the wavelet function and j and k are 
integers that scale and dilate the wavelet function. The 
factor „j" is known as the scale index, which indicates 
the wavelef's width. The location index k provides 
the position. The wavelet function is dilated by 
powers of two and is translated by the integer k. In 
terms of the wavelet coefficients, the wavelet equation 
is 

iU-l 

q?(x) = ^ gk^f20Qx - fc) 
k 

Where, gO, gl, g2 are high pass wavelet coefficients. 
Writing the scaling equation in terms of the scaling 
coefficients as given below, we get; 


JV-L 

0(a) = ^ hk J j20[2x — 
k 
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Figure 2: A 1-Dimensional DWT - Decomposition 
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Figure 3: One level wavelet decomposition of two 
dimensional data 



Figure 4: Sub bands in wavelet decomposition of two 
dimensional data 
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DWT is the multi resolution description of an image. 
The decoding can be processed sequentially from a 
low resolution to the higher resolution. DWT splits 
the signal into high and low frequency parts. The high 
frequency part contains information about the edge 
components, while the low frequency part is split 
again into high and low frequency parts. The high 
frequency components are usually used for 
watermarking since the human eye is less sensitive to 
changes in edges. In two dimensional applications, for 
each level of decomposition, we first perform the 
DWT in the vertical direction, followed by the DWT 
in the horizontal direction. After the first level of 
decomposition, there are 4 sub-bands: LL1, LH1, 
HL1, and HH1. For each successive level of 
decomposition, the LL Sub-band of the previous level 
is used as the input. To perform second level 
decomposition, the DWT is applied to LL1 band 
which decomposes the LL1 band into the four sub¬ 
bands LL2, LH2, HL2, and HH2. 

To perform third level decomposition, the DWT is 
applied to LL2 band which decompose this band into 
the four sub-bands - LL3, LH3, HL3, HH3. This 
results in 10 sub-bands per component. LH1, HL1, 
and HH1 contain the highest frequency bands present 
in the image tile, while LL3 contains the lowest 
frequency band. 

DWT is currently used in a wide variety of signal 
processing applications, such as in audio and video 
compression, removal of noise in audio, and the 
simulation of wireless antenna distribution. 
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Figure 5: 31evel discrete wavelet decomposition 

Wavelets have their energy concentrated in time and 
are well suited for the analysis of transient, time- 
varying signals. Since most of the real life signals 
encountered are time varying in nature, the Wavelet 
Transform suits many applications very well. As 
mentioned earlier, the wavelet equation produces 
different wavelet families like Daubechies, Haar, 
Coiflets, etc. The filter lengths and the number of 
vanishing moments for four different wavelet families 
are tabulated in Table 1. 


Wavelet Family 

Filters length 

Number of vanishing moments. N 

Haar 

2 

1 

Daubechies M 

1M 

M 

Coiflets M 

6M 

2M-1 

Symlets 

2M 

M 


Table 1: Wavelet families and their properties 


The term wavelet thresholding is explained as 
decomposition of the data or the image into wavelet 
coefficients, comparing the detail coefficients with a 
given threshold value, and shrinking these coefficients 
close to zero to take away the effect of noise in the 
data. The image is reconstructed from the modified 
coefficients. This process is also known as the inverse 
discrete wavelet transform. During thresholding, a 
wavelet coefficient is compared with a given 
threshold and is set to zero if its magnitude is less 
than the threshold; otherwise, it is retained or 
modified depending on the threshold rule. 
Thresholding distinguishes between the coefficients 
due to noise and the ones consisting of important 
signal information. 

Wavelet Denoising 

There are many techniques for image denoising which 
are using wavelet & based on wavelet shrinkage & 
wavelet thresholding. The computational advantage of 
using such estimation approaches is provided by 
algorithms of fast implementation. The smaller 
coefficients are instead eliminated, hence sparsifying 
the wavelet expansion. The simplest model of Image 
Degradation/Restoration Process: 



Noise 


i!(*j) 

Figure 6: Model of the image degradation/restoration 

process 

Wavelets provide a good expansion basis for the 
unknown function, because it can satisfy both 
parsimony & sparsity properties simultaneously; 
while the latter comes front the inherent distribution 
of the "energy" of the function over the coefficient 
vector, hence naturally selecting a subset of them as 
the most significant. This thesis work involves the 
denoising of intensity images with formats like JPG. 
Bto.lP. PNG & TIF[7]. 
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One of the most important properties of wavelets is 
that they provide unconditional bases for many other 
different signal classes. Hence, most of the signal 
information in wavelet expansions is conveyed by a 
relatively small number of large coefficients. The 
process of denoising using wavelets has been 
described by [14] & called “Wavelet Shrinkage”. In 
many systems, the concept of additive white Gaussian 
noise (AWGN) is used. This simply means a noise, 
which has a Gaussian probability density function & 
white power spectral density function (noise 
distributed over the entire frequency spectrum) & is 
linearly added to whatever signal used for analyzing. 

This work focuses primarily on thresholding as a 
means of denoising in the wavelet-domain due to 
some important characteristics present in wavelet 
coefficients. 

The noisy input signal can be thought of as the sum of 
the desired signal component (or true signal) & the 
additive white Gaussian noise (AWGN) with variance 
of a 2 & is shown in equation below. 

S(x,y ) = f(x,y ) + a 2 n(x,y ) 

Where; f ( x,y ) is desired signal & n (x,y) is noisy 
component. It has been shown that when the wavelet 
basis selected, matches well with the signal 
characteristics, a very few of the wavelet detail 
coefficients are influenced by the signal, while most 
of them are influenced by the noise. Therefore, an 
expression for the wavelet coefficients at each 
decomposition level can be described by the equation 
given as: 

y j =w j (Q + n j (Q .( 3 . 3 ) 

In addition, the desired signal coefficients are 
expected to be of larger magnitude when the SNR is 
not too small. Therefore, denoising is accomplished 
by thresholding wavelet coefficients, thereby 
eliminating noise-only coefficients & keeping the 
desired signal coefficients for reconstruction. 

III. PROPOSED METHODOLOGY 

In this part of the paper a detailed stepwise procedure 
is mentioned which is followed while denoising 
images using Orthonormal wavelets basis with eight 
vanishing moments (sym8) over four decomposition 
stages. The wavelet coefficients that lie on the same 
dyadic tree are well known to be large together in the 
neighborhood of image discontinuities. What can, 


thus, be predicted with reasonably good accuracy are 
the position of large wavelet coefficients out of 
parents at lower resolutions. However, getting the 
actual values of the finer resolution scale coefficients 
seem somewhat out of reach. This suggests that the 
best we can get out of between-scale correlations is 
segmentation between regions of large and small 
coefficients. In critically sampled orthonormal 
wavelet decomposition, the parent sub band is half the 
size of the child sub band. The usual way of putting 
the two sub bands in correspondence is simply to 
expand the parent by a factor two. Unfortunately, this 
approach does not take into account the potential non¬ 
integer shift caused by the filters of the DWT. We, 
thus, propose a more sophisticated solution, which 
addresses this issue and ensures the alignment of 
image features between the child and its parent. This 
idea comes from the following observation: Let 
LHj and LLj be, respectively, band pass and low pass 
outputs at iteration of the filter bank. Then, if the 
group delay between the bandpass and the low pass 
filters are equal, no shift between the features of Let 
LHj and LLj will occur. Of course, depending on the 
amplitude response of the filters, some features may 
be attenuated, blurred, or enhanced, but their location 
will remain unchanged. When the group delays differ, 
which is the general case, we, thus, propose to filter 
the low pass sub band LLj in order to compensate for 
the group delay difference with LHj. LLj is filtered in 
the three bandpass “directions” by adequately 
designed filters W HL , W HH & W LH providing aligned, 
i.e., group delay compensated, sub bands with 
HLj, HHj & LHj. Wavelet transform is appropriate 
for the digital image denoising methods. Because it 
can be transforming the images into sub-bands, each 
one contains certain frequencies that analyzed and 
processed easily [3], 

For better & easy understanding, a complete flowchart 
of the discussed methodology has been shown at the 
end of this chapter. The proposed algorithm steps are 
as follows: 


@ IJTSRD | Available Online @ www.ijtsrd.com | Volume - 2 | Issue - 6 | Sep-Oct 2018 


Page: 223 




International Journal of Trend in Scientific Research and Development (IJTSRD) ISSN: 2456-6470 



Figure 7: Proposed Denoising Algorithm 


IV. SIMULATION RESULTS & 

DISCUSSIONS 

PSNR is calculated for all the standard images with 
their noisy and denoised image counterparts, 
respectively. Hence I get good amount of comparison 
between the noisy and denoised image keeping the set 
standard image intact. 


PSNR stands for the peak signal to noise ratio. It is a 
term used to calculate the ratio of the maximum 
power of a test signal and the power of noise 
corrupted version of the test signal. Since most of the 
signals have a large dynamic range, PSNR is 
generally represented in terms of the logarithmic 
decibel (dB) scale. It is most commonly used as a 
measure of quality of reconstruction in image 
compression etc. It is calculated as the following: 


^—1771 — 1 i7l— 1 

MSE = -) ) 11 /( 6 ;) - K(i,j 

mn 4_i i= o Zj; =0 


j) IP 


(MSEf\ 

PSNR = 10., <o Sl0 ( w ) 

Where I and K are the original and noisy / denoised 
image, respectively. MAXj is the maximum pixel 
value of the image under test. For an image having 8 
bits per sample, pixels representation, this is 
equivalent to 2 8 = 255. 

In simulation we have done the simulations with 
Orthonormal wavelets basis with eight vanishing 
moments (sym8) over four decomposition stages. 
Standard test image ‘Boat’ is taken for verification of 
our proposed algorithm. We have performed 
simulations for AWGN for a range from 5 dB to 30 
dB. 


Original Noisy: PSNR = 34.15 dB Denoised: PSNR = 36.73 dB Mean SSIM: 0.9829 



100 200 300 400 500 100 200 300 400 500 100 200 300 400 500 

Figure 8: Test image ‘Boat proposed method using 
‘Coiflet’ wavelet at 5 dB AWGN noise 


Original Noisy: PSNR = 28.13 dB Denoised: PSNR = 32.94 dB Mean SSIM: 0.9558 



Figure 9: Test image ‘Boat’ proposed method using 
‘Coiflet’ wavelet at 10 dB AWGN noise 


Original Noisy: PSNR = 24.61 dB Denoised: PSNR = 30.90 dB Mean SSIM: 0.9262 



100 200 300 400 500 100 200 300 400 500 100 200 300 400 500 

Figure 10: Test image ‘Boat proposed method using 
‘Coiflet’ wavelet at 15 dB AWGN noise 


Original Noisy: PSNR = 22.11 dB Denoised: PSNR = 29.51 dB Mean SSIM: 0.8972 



100 200 300 400 500 100 200 300 400 500 100 200 300 400 500 

Figure 11: Test image ‘Boat’ proposed method using 
‘Coiflet’ wavelet at 20 dB AWGN noise 





Original Noisy: PSNR = 20.17 dB Denoised: PSNR = 28.48 dB 


Figure 12: Test image ‘Boat’ proposed method using 
‘Coiflet’ wavelet at 25 dB AWGN noise 


Mean SSIM: 0.8697 
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Original Noisy: PSNR = 18.59 dB Denoised: PSNR = 27.67 dB Mean SSIM: 0.8442 


100 200 300 400 500 100 200 300 400 500 100 200 300 400 500 

Figure 13: Test image ‘Boat’ proposed method using 
‘Coiflet’ wavelet at 30 dB AWGN noise 


IV.2 SIMULATION RESULTS SUMMARY 


In simulation we have taken Symlet, Coifilet etc 
orthogonal wavelets. The summary of experiments 
using proposed method is listed in table 2. 


Noise 

Standard 

Deviation 

Noisy Image Parameters 

Denoised using Proposed OWT- SURELET 
by Coiflet Wavelet 

<r (dB) 

Noisy image 
PSNR (dB) 

Noisy image 
MSE 

PSNR (dB) 

MSE 

SSIM 

5 

34.15 

25.00 

36.73 

13.80 

0.9829 

10 

28.13 

100.00 

32.94 

33.02 

0.9558 

15 

24.61 

225.00 

30.90 

52.86 

0.9262 

20 

22.11 

400.00 

29.51 

72.73 

0.8972 

25 

20.17 

625.00 

28.48 

92.24 

0.8697 

30 

18.59 

900.00 

27.67 

111.21 

0.8442 


Table 2: Test Image ‘Boat’ Simulation Results 


PSNR Comparison 
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Figure 14: PSNR comparison for test image ‘Boat’ 

MSE Comparison 
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Figure 15: MSE comparison for test image ‘Boat’ 

SSIM using Proposed Methodby Coiflet Wavelet 
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Figure 16: SSIM using Proposed method by Coiflet 
Wavelet for test image ‘Boat’ 





IV.3 SIMULATION RESULTS COMPARISON 


In Image Denoising Using Discrete Wavelet 
Transform and Edge Information [1], denoising with 
two different methods has been performed for various 
test images. A comparison for standard ‘Boat’ image 
is listed in table 3. 



Table 3: Test Image Boat Simulation Results 
Comparison 

PSNR Comparison 


35 



Noise Standard Deviation a (dB) 

■ Noisy image ■ Spatial Gaussian filter 5x5 

■ Conventional filtering using DWT and soft threshold ■ Denoising with one threshold for edge detection [ 1] 

■ Denoising with two thresholds for edge detection [1] ■ Proposed method with OWT SURELET 

Figure 17: PSNR comparison graph for test image 
‘Boat’ 

V. CONCLUSION 

In this paper image denoising techniques for the 
AWGN corrupted has been given. A new algorithm 
based on the orthonormal wavelet transform (OWT) is 
developed. In this work images corrupted by AWGN 
are denoised. Simulation results shows that proposed 
method using Orthonormal wavelets for different 
values of noise Standard Deviation a in dB 
outperforms other available methods. Also Coiflet 
Wavelet performs better than Symlet, Haar & 
Daubechies wavelets. The proposed Orthonormal 
wavelet transform (OWT) method has minimum 
Mean Square & highest PSNR with Coiflet wavelets. 
Simulation results shows that denoised image is 98.29 
% similar for 5 dB noise standard deviation and 
84.42% similar for 30 dB noise standard deviation. 
The proposed method has 1.35 dB & 4% 

improvement for a =10 dB, 2.08 dB & 7% 
improvement for a =20 dB & 2.26 dB & 9% 
improvement for a =30 dB as compared to denoising 
with two thresholds for edge detection [1], 
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